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SUMMARY 

An analysis vas made to determine the maximum power plant thrust-to-weight 
ratio possible with a single-cavity vortex gaseous reactor in which all the hy¬ 
drogen propellant must diffuse through a fuel-rich region. An assumed radial 
temperature profile was used to represent conduction, convection, and radiation 
heat-transfer effects. The effect of hydrogen property changes due to dissocia¬ 
tion and ionization was taken into account in a hydrodynamic computer program. 

It is shown that, even for extremely optimistic assumptions of reactor crit¬ 
icality and operating conditions, such a system is limited to reactor thrust-to- 
weight ratios of about 1.2x10-3 for laminar flow. For turbulent flow, the maxi¬ 
mum thrust-to-weight ratio is less than 10 *" 3 . These low thrusts result from the 
fact that the hydrogen flow rate is limited by the diffusion process. The per¬ 
formance of a gas-core system with a specific impulse of 3000 seconds and a pow- 
erplant thrust-to-weight ratio of 10“ 2 is shown to be equivalent to that of a 
1000-second advanced solid-core system. It is therefore concluded that a single¬ 
cavity vortex gaseous reactor in which all the hydrogen must diffuse through the 
nuclear fuel is a low-thrust device and offers no improvement over a solid-core 
nuclear-rocket engine. 

To achieve higher thrust, additional hydrogen flow must be introduced in 
such a manner that it will bypass the nuclear fuel. Obviously, such flow must be 
heated by thermal radiation. An illustrative model of a single-cavity vortex 
system employing supplementary flow of hydrogen through the core region is 
briefly examined. Such a system appears capable of thrust-to-weight ratios of 
approximately 1 to 10. For a high-impulse engine, this capability would be a 
considerable improvement over solid-core performance. Limits imposed by thermal 
radiation heat transfer to cavity walls are acknowledged but not evaluated. Al¬ 
ternate vortex concepts that employ many parallel vortices to achieve higher hy¬ 
drogen flow rates offer the possibility of sufficiently high thrust-to-weight 
ratios, if they are not limited by short thermal-radiation path lengths. 


INTRODUCTION 


The use of a reactor employing a fissionable fuel in a gaseous state has re¬ 
ceived considerable attention as a propulsion unit for a nuclear rocket. Such a 
device offers the attractive possibility of improved performance by increasing 
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the maxim-urn propellant-exhaust temperature "beyond the limit imposed by the melt¬ 
ing point of solid-fuel elements in a heat-transfer nuclear rocket. An elemen¬ 
tary consideration of hydrogen-propellant flow rates required for reasonable ac¬ 
celerations (e.g., 10 g*s or greater) and nuclear-fuel-to-hydrogen-mass ratios 
required for criticality leads immediately to the conclusion that the fuel con¬ 
sumption of a "straight-through" gaseous reactor would be prohibitive. To cir¬ 
cumvent this problem, certain flow schemes have been considered that employ hy¬ 
drodynamic forces to increase the holdup time of the heavy fuel molecules rela¬ 
tive to that of the light hydrogen molecules as the mixture passes through the 
reactor cavity. This increased holdup time decreases the required mass—flow rate 
of fuel for a given hydrogen flow rate. 

Since 1958, the principal effort in gaseous-reactor research has been di¬ 
rected toward a study of a vortex system. In such a device, the propellant (hy¬ 
drogen) and fissionable fuel (both uranium 255 and plutonium 239 have been con¬ 
sidered) are injected at the outer wall of a reactor cavity with a velocity com¬ 
ponent in the tangential direction. The gases swirl radially inward and are fi¬ 
nally exhausted axially. The pressure gradient attendant to this flow field in¬ 
duces pressure diffusion that tends to move the heavier component (fuel) toward 
the outer wall relative to the light one (hydrogen). This effect is modified by 
concentration diffusion and radial inflow forces. The steady-state result is 
that an annular "cloud" forms that is rich in heavy gas. The hydrogen gas dif¬ 
fuses through this region and acquires the nuclear energy release directly by 
molecular collisions. This general flow pattern is illustrated In figure 1. 

Predictably, the emphasis of gaseous-vortex-reactor research has been on 
the hydrodynamics of the system (refs. 1 to 3). The obvious reason is that if 
the flow pattern does not result in a critical fuel concentration for acceptable 
hydrogen-to-fuel flow-rate ratios further considerations are not necessary. 

In addition to the basic separation process, other limits of performance 
capability exist. Some of these performance limits are inherent in any exter¬ 
nally moderated gaseous-reactor system, and others apply specifically to a vortex 
concept. Certain general limitations have been determined for the spectrum of 
nuclear-rocket systems from all solid-core to all gas-core systems (refs. 4 to 6). 
A number of recent investigations (refs. 7 to 10) provide a basis for estimating 
the effects of turbulence in a vortex flow field. Based on certain simplifying 
assumptions, an analysis has been made of the performance of a turbulent gaseous- 
reactor system (ref. 11). Some effects of thermal radiation (refs. 12 to 14) and 
criticality characteristics (refs. 13 and 15) have also been reported, as well as 
overall system studies (ref. 16). 

In view of the combined scope of these various studies, it seems that an 
evaluation of the most optimistic performance capability of a single-cavity vor- 
tex gaseous nuclear rocket is now possible and of value. The goal here Is to es¬ 
tablish the maximum performance capability and not the most probable. Further, 
there are, in a sense, two maximum performance levels; one is associated with 
laminar flow, the other with turbulent flow in the vortex. While it is impossi¬ 
ble to say that some future technique will not result in a laminar high-strength 
vortex, published experimental information indicates that turbulence Is most 
likely to exist. 
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Finally, a -word about the criterion of maximum performance to he used is in 
order. The system considered herein is one in -which all of the nuclear fuel is 
a gas and is in the reactor cavity. This situation leads to a fixed maximum 
specific impulse for all gaseous-reactor systems and eliminates the need to con¬ 
sider mission requirements and overall vehicle performance. The hydrodynamics of 
the system are such that an increase of hydrogen flow rate results in a decrease 
of average fuel concentration in the reactor cavity; hut, to maintain critical¬ 
ity, a larger and heavier reactor is then required. To balance these two ef¬ 
fects, the criterion of performance used in this report is the thrust-to-weight 
ratio of the reactor. 

Hydrogen property variations are included in the hydrodynamic computer pro¬ 
gram for the IBM 7090 digital computer given in appendix B hy Muriel B. Elan. 


PERFORMANCE LIMITATIONS 


It is possible to establish certain limits on a gaseous-vortex-reactor sys¬ 
tem before making a detailed analysis of the fluid mechanics. Some of these 
limitations are general in nature and apply to any externally moderated reflected 
gaseous-reactor system; others apply only to the vortex concept. Some of the 
limits discussed are arbitrary; in these instances, opt imi stic values are chosen 
to establish the maximum bounds of performance. 


GENERAL LIMITATIONS 
Fuel Consumption 

The parameter governing fuel consumption is the hydrogen-to-fuel mass-flow- 
rate ratio. An upper limit will be established by hydrodynamic and nuclear con¬ 
siderations. If the ratio is too large, fuel concentrations necessary for criti¬ 
cality are not possible. A lower limit on this ratio Is somewhat less clearly 
defined, at least from economic considerations. A ratio of the order of 1000 to 
1 has been used as an acceptable value in past considerations. For the purpose 
of this report, a hydrogen-to-fuel mass-flow-rate ratio of 100 to 1 is taken ar¬ 
bitrarily. Values ranging from 50 to 1 up to 10, 000 to 1 were investigated to 
show the effect of this parameter. 


Propellant Molecular Weight 

Consumption of nuclear fuel is not the only factor Involved in the lower 
limit of propellant-to-fuel flow ratio. As the fuel-to-hydrogen atom ratio in¬ 
creases, the average molecular weight of the exhaust gas obviously increases. An 
upper limit on fuel concentration is reached when the increase In the average 
molecular weight begins to decrease the specific impulse significantly. The ex¬ 
haust gas average molecular weight M can be easily expressed In terms of the 
hydrogen-to-fuel mass-flow-rate ratio and the fuel-to-hydrogen molecular-weight 
ratio: 
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For the situations considered here, (w^/wg) (Mg/M-^) is sufficiently large so that 



The 100 to 1 mass-flow*-rate ratio selected from economic considerations will re¬ 
sult in approximately a 1-percent increase in propellant molecular weight. The 
resulting l/2-percent decrease in specific impulse certainly does not represent 
a serious limitation on performance, and the mass-flow-rate ratio of 100 to 1 
will give optimistic results. (All symbols are defined in appendix A.) 


Nucleonics 

There are two limitations that arise from the nucleonics of gaseous reac¬ 
tors. The first of these is the pressure required for criticality. With all 
other parameters fixed, virtually any reactor scheme can he made critical by sim¬ 
ply operating at whatever system pressure is required for criticality. The 
choice here, again arbitrary, will be a maximum pressure in the reactor cavity of 
500 atmospheres. Since the contemplated operating pressure of current solid-fuel 
reactors is considerably smaller, this choice must be considered optimistic. 

A more subtle l i mitation exists due to gamma and neutron heating of the 
solid materials surrounding the reactor cavity. This effect has been considered 
in detail elsewhere (ref. 4) and is only briefly reiterated here. The fraction 
of the total power generation, which results from gamma and neutron heating, and 
thermal radiation to the cavity walls must be removed from the solid regions by 
the hydrogen before it enters the reactor cavity. The resulting enthalpy rise of 
the hydrogen in the solid region of a gaseous reactor clearly cannot exceed the 
total rise possible in a solid-core-reactor system. This, then, fixes the maxi¬ 
mum possible enthalpy gain available to the hydrogen in a gas-core system, as 
explained in the following paragraph. 

It is assumed herein that the combined effects of nuclear and thermal radia¬ 
tion result in the deposition of 10 percent of the total power in the solid re¬ 
gion. Further, it is assumed that the hydrogen enters the reactor cavity at 
5000° B. These conditions result in a maximum possible exhaust temperature of 
15,000° to 20,000° B, depending on pressure, and a specific impulse of approxi¬ 
mately 3000 seconds. The only possibility of obtaining a higher specific impulse 
is through the use of a radiator to dispose of the additional solid-region heat 
(ref. 4). 
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Heat Transfer to Walls 


Two aspects of gaseous-reactor heat transfer are mentioned here, although no 
resulting limitations are included in this study. The first limit is imposed hy 
the allowable cavity wall temperature. Any system that proposes to contain gas 
at temperatures of the order of 15,000° R is faced with thermal-radiation prob¬ 
lems; these have been treated in other investigations (refs. 12 to 14). Another 
heat-transfer problem is associated with exhaust-nozzle cooling. Although rela¬ 
tively unexamined to date, this problem may well place a serious limitation on 
any gaseous-reactor system. 


VORTEX-SYSTEM LIMITATIONS 
Geometry 

The performance limitation imposed by reactor geometry arises from the fact 
that the separation process improves with decreasing mass flow of hydrogen per 
unit cavity length. The combination of criticality and hydrodynamic considera¬ 
tions suggests that the reactor should be composed of many small-diameter tubes 
in order to obtain a high hydrogen mass flow (refs. 1 and 17); however, the inlet 
headers, exhaust ducting, and short thermal radiation path lengths associated 
with small diameter, large length~to-diameter ratio (L/L) vortex tubes pose ex¬ 
tremely severe heat-transfer problems. All the findings of this report, as the 
title implies, are for a single-cavity vortex gaseous reactor. 


Hydrodynamics 

For given values of critical fuel density and maximum operating pressure, 
the hydrodynamic relations prescribe a maximum hydrogen flow rate per unit cavity 
length. The assumptions utilized in reference 1 gave a value of approximately 
0.01 pound per second per foot. The two flow regions to be considered herein, 
laminar and turbulent, will result in one maximum for each. This maximum hydro¬ 
gen flow rate, however, will not necessarily yield the highest thrust-to-weight 
ratio. Two other limitations imposed by hydrodynamics are: (l) axial choking 
at the vortex exit, and (2) nonequilibrium nozzle expansion. The limit of choked 
flow in the vortex exhaust is included in the analysis of this report, and equi¬ 
librium expansion is assumed. 


ANALYSIS 


An analysis of the hydrodynamics, vortex axial choking, criticality, and 
performance characteristics of a vortex gaseous nuclear rocket is presented 
herein, in that order. The general development of the hydrodynamic relations 
follows that presented in reference 2. Hydrogen properties and the hydrodynamic 
computer program are discussed in detail in appendixes B and C. 


5 



ASSUMPTIONS 


The assumptions utilized in the following analysis are listed here for con¬ 
venience. The more important ones are also noted at the point in the analysis 
where they are introduced. 

(1) b/b& « b/bz s= 0 (except in vortex core where axial mass flow varies 

linearly with axial position) 

(2) Steady-state conditions exist 

(3) Perfect gas law is obeyed 

(4) Hydrogen gas species are present in equilibrium amounts 

(5) External body forces are negligible 

(6) Thermal diffusion is small compared with pressure and concentration 
diffusion 

(7) Radial velocities are small compared with tangential velocities 

(8) p and pe do not vary with radial position 

(9) Laminar diffusion coefficient is given by Gilliland expression 

(10) Turbulent flow is described by laminar equations if j 1 is replaced by 
p€ m0 and (B-L 2 + e mR ) is substituted for 

(11) Turbulence affects concentration diffusion but not pressure diffusion 

(12) Combined effect of molecular collision heating, conduction, and radia¬ 
tion heat transfer is adequately represented by an assumed radial tem¬ 
perature profile t(x) 

(13) Reactor criticality is unaffected by a nonuniform fuel distribution or 
presence of structural material 

(14) Equilibrium expansion occurs in exhaust nozzle 

(15) Ten percent of total power generation is deposited in moderator- 
reflector region 

(16) The maximum possible hydrogen peripheral Mach number is 1.0 

(17) Fuel gas does not ionize 

(18) Axial choking occurs at density conditions in vortex at exit radius 

(19) Maximum permissible pressure at vortex perimeter is 500 atmospheres 

(20) Temperature of hydrogen entering reactor cavity is 5000° R 
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HYDRODYNAMICS 


Laminar Flow 

The Navier-Stokes equations of motion for compressible vortex flow are as 
follows: 



These equations are consistent with assumptions (l), (2), (7), and (8). 

In lieu of an energy equation, the static-temperature profile will he 
represented by an assumed distribution (assumption (12)) of the form 

t(x) = a 1 + a 2 x + a 3 x 2 + a 4 x 5 + a 5 x 4 (3) 

The constants are chosen such that temperature increases monotonically 
from 5000° R at the outer radius to a maximum value at the center of the 
vortex and that the temperature gradient is essentially zero at the center. 

In addition, the temperature increases to near its maximum value at a radius 
ratio larger than that of the maximum uranium concentration for all cases con¬ 
sidered. This would actually occur because conduction and radiation would 
tend to increase the hydrogen temperature before it entered the fuel-rich 
region. 

An alternate approach would be to assume that the volumetric heat- 
generation rate is proportional to the concentration of fissionable gas 
(ref. l). Since a rigorous treatment of radiative heat transfer would intro¬ 
duce considerable complexity to an already cumbersome computer program, the 
compromise represented by equation (3) is considered entirely adequate for 
this analysis. 

The continuity equation for the annular region r Q ^ r ;> r n is 

^ Pur = 0 (4) 

The assumption of a perfect gas gives 


- PRt 
~ M 


(5) 


as the equation of state. 

By virtue of assumptions (5) and (6), the diffusion equation is 
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a l n 2 


+ n i n 2^ M 2 ~ ^ 1 d£ 
dr np p dr 


where 


n 2 - 


( 6 ) 


A laminar Reynolds number is defined in terms of the radius and radial 
velocity as 


Re- 


pur 


(?) 


The negative sign is used so that a, positive Reynolds number is associated 
with flow that is radially inward. The continuity equation and assumption ( 8 ) 
result in a Reynolds number that does not vary with radial position. 

Integration of the continuity equation for the annular region yields 


-pur = 2 ^ = 


2 * 


( 8 ) 


where w is the total mass flow per unit axial length. In the core of vortex 
r n > r > 0 ^ the assumption of constant mass flow per unit area in the axial 
direction gives the continuity relation (ref. 10 ) as. 


pur = P n %r n (i-) 


O) 


The tangential velocity profile is obtained by integrating equation (2) 
with the boundary conditions v(r Q ) = v Q and v( 0 ) = 6 ; 


-(*)v 


where x is the radius ratio r/r Q , and f is a function of 


Re 7 


In the vortex core, \j r is given by 


♦ M 


and 


r 

x n 

-Ee z /2(x/x n ) 2 

1 - e 



I Ri 72 

r 


1 - e 




j 


( 10 ) 


( 11 ) 


% = [ + H^7)( e Re z?_ J( xRn 2 -1 


(lla) 


In the annular region, \jr is given by 
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( 12 ) 


■where 





(i - c ) 2 





(12a) 


and X n / Y Q is given "by equation (lla). 


It is of interest to note that -when j i 0, the function \|r -► 1, and equa¬ 
tion (10) gives the usual l/r variation of tangential velocity for a potential 
vortex. In fact, the l/r profile is closely approximated for any Re^ greater 

than about 10. 

The pressure gradient from equation (l) may he -written in the form 


_a/p\ »i M i + y2M 2 K 2 /p\ 

dx\p*/ gRtx \P*/ 


The pressure is normalized to p*, its value at some radius ratio x*. In 
this analysis, x* is the radius ratio at ■which the mole fraction of the heavy 
component y 2 is at a maximum yg. The numerical integration procedure will 

begin at x* and proceed to x = 1, then go from x* to x n , and, finally, 
from x n to 0. The computer program is discussed in detail in appendix B. 


The diffusion equation may be written in terms of y 2 using equations (l), 
(5), (6), and (8) as follows: 


where 


dy 2 a 2 (by 2 - l) . P(y 2 - y|) 
dx ~ x x 3 


w 2 

a2 2*M 2 D 12 n 


(Mg - M 1 )Vq11t 

gRt 


b = 1 + 


^2 

W 2 M x 
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The product of the laminar diffusion coefficient and the molecular density 
is calculated from the Gilliland relation: 

l/2 

2 . 62 xlO- 6 -Vt(j|- + j|) 

D 12 n = ( v l/3 + v l/3)2 

The molecular volumes of hydrogen and uranium were taken to he 14. 3 and 105 cubic 
centimeters per gram mole, respectively. The other hydrogen properties, molec¬ 
ular weight, specific heat, ratio of specific heats, and composition were allowed 
to vary with vortex radius as a function of temperature and pressure. The hydro¬ 
gen property variation is discussed in appendix C. 

The mole fraction of the heavy component y 2 is related to mass-flow and 
radial-velocity ratios in the following manner: 


y 2 « 


1 


1 + 


■Wi M2 U2 
v 2 M x u-l 


(15) 


The system of equations required to describe the hydrodynamic separation 
process is now complete and is given by equations (3), (5), either (8) or (9), 
(10), (13), and (14). The two nonlinear differential equations were solved by 
the Runge-Kutta numerical method. The boundary condition for equation (13) is 
p - p* at x = x*. The boundary condition for equation (14) is dy 2 /dx* = 0. 

The implications of this boundary condition are discussed in RESULTS AND DISCUS¬ 
SION. 


Total temperature and pressure are calculated from 

2 


T - t + 


2c p gJ 


(16) 


P = p 


(r - ^ 

(yA + y 2 M 2 )v2 l 

V 2 ) 

rgRt J 


r/(r-i) 


(17) 


Turbulent Flow 


The hydrodynamic relations for a turbulent vortex are obtained by substitut¬ 
ing pe mo for |-i and D 12 + for D 12 into the laminar equations with one 

exception in the diffusion expression. For turbulent flow, the Reynolds number 
is defined as 
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( 18 ) 


The tangential velocity profile is still given by equation (10), hut i|r is now 
evaluated with a turbulent Reynolds number. The eddy diffusivity for momentum 
transfer e mo can be calculated from equation (18) for a given value of Re^. 

For all the turbulent-flow considerations in this report, e mQ /= 1. 

It will be assumed herein, as it was in reference 2, that equation (14) is 
valid for turbulent diffusion if (D-j_ 2 + e ma ) is substituted for 1 12 in the 
concentration-gradient term of equation (6). The coefficient of the pressure- 
diffusion term is unchanged, but that does not imply that pressure diffusion is 
unaffected by turbulence. The coefficient is unchanged, but the effect of tur¬ 
bulence is exhibited in the pressure gradient. 

With this assumption, the turbulent-diffusion equation becomes 


ay 2 Dig f-ogCbyg 1) + p(y 2 - y|) 

^ d 12 + € ma x x 3 


(19) 


The functions Og, b, and (3 are the same as those given for equation (14). 


AXIAL CHOKING 

An upper limit on hydrogen flow rate is placed by hydrodynamic choking at 
the vortex axial exit. The hydrogen flow rate per unit length required for 
choked flow w ch is given by 


v s A nPn 
v ch L~ 


( 20 ) 


where v g , the velocity of sound in hydrogen, and the density p n are evaluated 
at the vortex-exit radius r n . This approximation will result in an optimistic 
(high) value for w ch ; the average density in the vortex exit will be less than 
the value at r n because of the pressure gradient in the core. 

For temperatures and pressures in this study, y for an equilibrium mixture 
of hydrogen species varies less than 1 percent from 1.05. With this value and 
A n = 2*r n , equation (20) can be written in the following form for t n = 19, 860 R: 

where Mn M is the average molecular weight of the propellant at p n and T n . 
j., n 
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CRITICALITY 


The critical mass values used In this report are taken from the results of 
a two-dimensional analysis of gas-core cavity reactors (ref. 15). Critical den¬ 
sity as a function of reactor radius for cavity length-to-diameter ratios of 
1/2, 1, and 2 was computed herein from cylindrical reactor critical mass esti¬ 
mates “based on a “buckling analogy presented in reference 15. These criticality 
requirements are for a reactor cavity uniformly filled with uranium 235 and com¬ 
pletely surrounded by a 70° F heavy water (D 2 0) moderator-reflector region. 

This choice of moderator material admittedly results in a criticality re¬ 
quirement that is more optimistic than realistic. Estimates of critical density 
for two other reactor systems were obtained by using multiplying factors obtained 
from reference 15. One of these systems is a uranium 235 fueled cavity sur¬ 
rounded by graphite at 5300° F, and the other is a plutonium 239 fueled cavity 
with a 5300° F graphite moderator. 


PERFORMANCE 

Thrust-to-reactor weight is used as the criterion of reactor performance. 
This ratio Is given by: 


JL _ v i l:i: sp 

W R W R 


( 22 ) 


It is assumed that the reactor weight is that of the moderator reflector. In 
equation (22), L is the length of the reactor cavity; it does not include the 
end reflectors. 


RESULTS AND DISCUSSION 
GENERAL PERFORMANCE CHARACTERISTICS 
Comparison With Solid-Core Performance 

It is worthwhile to investigate some general characteristics of gas-core 
nuclear rockets that are independent of the hydrodynamics of a specific system. 
The trade-off between specific Impulse and powerplant thrust-to-weight ratio, for 
example, is of interest because a gas-core system may provide a higher impulse 
but lower thrust-to-welght ratio than a solid-core heat-transfer rocket. 

Figure 2, which is taken from information in reference 18, shows the effect 
of powerplant thrust-to-weight ratio on the specific impulse required for a 
seven-man 460-day Venus mission. The parameter here Is the ratio of vehicle 
weight for a system under consideration to that of an advanced (l gp = 1000 sec) 

solid-core system. The curves of figure 2 show that a gaseous-reactor scheme 
that provides a specific impulse of 3000 seconds must have a powerplant thrust- 
to-weight ratio of about 10“ 2 in order to perform the same mission as a solid- 
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core rocket of the same gross weight. A thrust-to-weight ratio of approximately 
2 with a 3000-second impulse would allow the vehicle weight to he one-half that 
of a solid-core system. So, the general conclusion here is that a gas-core sys¬ 
tem must offer reactor thrust-to-weight ratios greater than 10“ 2 in order to 
merit consideration. 

The performance advantage of a 3000-second-impulse engine with a thrust-to- 
weight ratio greater than 10"^ does not have to he utilized in terms of decreased 
vehicle weight. Figure 3 shows the decrease in trip time that is possible be¬ 
cause of increased specific impulse for a constant (1.7x10® lb) vehicle weight. 
The 460-day trip time required by a 1000-second solid-core rocket can be reduced 
significantly to 150 days for a 3000-second gas-core system. This gain over a 
solid-core rocket is probably of more value than a gain in terms of decreased 
vehicle weight. 


Attainable Specific Impulse 

Because a gaseous nuclear rocket can potentially produce extremely high 
propellant temperatures, it is of interest to consider briefly what specific 
impulses are possible with such systems, though a fixed impulse of 3000 seconds 
is used in this study. Since specific impulse is simply a measure of energy con¬ 
tent, specification of the reactor temperature and pressure determines the equi¬ 
librium impulse available from hydrogen exhausting to a vacuum. 

Figure 4 shows the equilibrium specific impulse in vacuum for temperatures 
from 5000° to 35,000° R and reactor pressures of 1, 10, and 100 atmospheres, as 
calculated from enthalpy values given in reference 19. As pointed out in PER¬ 
FORMANCE LIMITATIONS, the fraction of the total enthalpy rise that the hydrogen 
must acquire to remove the energy deposited in the solid regions of a reactor 
limits the total enthalpy rise possible. This parameter cp is shown in figure 4 
for an assumed hydrogen temperature of 5000° R at the inlet to the reactor cav¬ 
ity. If 10 percent of the total reactor power must be removed from the solid 
regions, the maximum hydrogen exhaust temperature is between 17,000° and 
19,000° R. This gives a specific-impulse range from 2900 to 3200 seconds. For 
a cp of 20 percent, hydrogen exhaust temperatures are limited to 7500° to 
10,000° R, with impulses from 2050 to 2275 seconds. Though values of 15 to 20 
percent solid-region heat generation are probably more realistic, a value of 
10 percent is assumed possible, and a constant specific impulse of 3000 seconds 
was used herein for all performance calculations. 

A vortex gaseous reactor enjoys one advantage as a result of the large ra¬ 
dial pressure gradient attendant to a vortex flow field. The radial decrease in 
pressure allows dissociation and ionization of the hydrogen to occur at a lower 
temperature (fig. 5). For a given wall pressure, the average exhaust tempera¬ 
ture necessary to obtain an enthalpy rise that is 10 times that of a solid-core 
rocket (operating at the same wall pressure) decreases at lower exhaust-nozzle 
pressure. For example, a uniform reactor pressure of 100 atmospheres would re¬ 
quire a hydrogen average exhaust temperature of 16, 750° R to give a specific im¬ 
pulse of around 3000 seconds. A vortex reactor operating at a wall pressure of 
100 atmospheres and an exhaust pressure of 0.1 atmosphere would only require an 
exhaust temperature of 14,950° R for the same impulse. 
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Reactor Criticality 


The criticality requirements of an externally moderated cavity reactor are 
shewn in figure 6. These curves were computed from results presented in refer¬ 
ence 15 and are for cylindrical reactors with side and end moderator reflectors 
100-centimeters thick. The three solid-line curves of critical density against 
cavity radius are for a uranium 235 fuel and a 70° F DgO moderator. For high- 

temperature operation, heavy water is obviously an unrealistic choice hut was 
used for the calculations of this report in order to assure an optimistic evalua¬ 
tion of performance. The criticality requirements of more practical systems 
that utilize 5300° F graphite as a moderator and either plutonium 239 or uranium 
235 as a nuclear fuel are shown for comparison, hut were not used. The curves in 
figure 6 are for a reactor cavity uniformly filled with fuel, and it is assumed 
here that no increase in critical mass would result from a nonuniform distribu¬ 
tion, though in fact this would he true. 


Reactor Geometry 

It is also possible to determine some general effects of reactor geometry 
on rocket performance that are independent of a detailed hydrodynamic analysis. 
Figure 7 shows reactor weight (assumed to he only that of the moderator) as a 
function of reactor-cavity radius and length-to-diameter ratio. The upper limits 
of these curves result from the fact that cavity radii greater than 300 centi¬ 
meters and/or reactor weights greater than 10 6 pounds are not considered herein. 

For a fixed hydrogen flow rate per foot of vortex length, an increase in 
reactor length results in an increase of both thrust and reactor weight. The 
solid curves in figure 8 show this geometry effect for cavity length-to-diameter 
ratios of l/2, 1, and 2. For a given cavity radius, a length-to-diameter ratio 
of 2 gives the largest thrust-to-weight ratio. It will he shown subsequently 
that values of length-to-diameter ratio greater than 2 yield only marginal im¬ 
provement of performance. 

There is, however, another factor to consider here. Since geometry also af¬ 
fects the criticality requirement, the reactor radius and the length-to-diameter 
ratio cannot he varied independently for a given average fuel density. The 
dashed lines in figure 8 show how vortex radius and length-to-diameter ratio must 
he varied in order to maintain a constant critical density. The conclusion to he 
drawn from the curves of figure 8 is that the more favorable performance will re¬ 
sult from a length-to-diameter ratio of 2 and the smallest radius necessary to 
maintain criticality. 


VORTEX HXDRODYEAMIC CHARACTERISTICS 
Parameters Considered 

For a given set of input flow and geometry parameters the hydrodynamic equa¬ 
tions give the radial distribution of the fissionable gas throughout the vortex. 
This distribution is then used to obtain: (l) average fuel density, (2) reactor 
size (for a given L/D), and (3) reactor weight. A description of the computer 
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program, format is given in appendix B. The input parameters investigated are as 
follows * hydrogen flow rate, uranium flow rate, peripheral tangential velocity, 
vortex-exit-radius ratio, radius ratio of maximum uranium mole fraction x , ra¬ 
dial Reynolds number (12 for laminar flow, 2.5 for turbulent flow), and radial 
temperature profile. The input static pressure at x* is not a primary param¬ 
eter; it was adjusted in every case to give a pressure at the outer radius of 
500 atmospheres. 


Effect of Temperature Profile 

Two assumed temperature profiles were investigated and are shown in figure 
9(a) Both start at 5000° R and increase with decreasing radius to a maximum 
value of 15,000° R, and of 20,000° R at the center of the vortex. The 20,000 R 
profiles was used in some of the initial parameter variations, hut the best 
performances, both turbulent and laminar, were obtained with the 15,000 R pro¬ 
file. 


The effect of the temperature profile on the uranium distribution for a 
typical laminar flow situation is shown in figure 9(b). The effect is ^slight 
and in the expected direction. The lower temperature results in a maximum con¬ 
centration that is approximately 20/15 higher. The increase of the average con¬ 
centration in the entire vortex is, of course, considerably less than that of 
the maximum. 


Effect of Hydrogen Plow Rate 


i 
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Laminar flow. - The effect of hydrogen flow rate on the separation process 
is shown in figure 10 for laminar flow. The ordinate, the ratio of uranium to 
hydrogen density, is a measure of how much "separation” has occurred. Ho separa¬ 
tion implies that the uranium streamlines through the vortex are the same as 
those for hydrogen, and so the density ratio of uranium to hydrogen is constant. 
This situation occurs for relatively high flow rates. As the hydrogen flow ra e 
is decreased, the amount of separation increases. 


Figure 10 illustrates the basic limitation of a single-cavity vortex system 
where all of the hydrogen must diffuse through the fuel region. The diffusion 
rates limit the hydrogen flow rates to small values for conditions under which 
separation will occur. The curves in figure 10 also show why the maximum hydro¬ 
gen flow rate for a given average uranium density does not necessarily give the 
maximum thrust per reactor weight. High flow rates give high thrust butJhe 
reactor critical size, and therefore the weight, becomes infinite. For extremely 
small flow rates, the reactor weight is reasonable, but the thrust is low. This 
effect is shown in figure 11 for laminar flow and a re act or-cavity length-to- 
diameter ratio of 2. There is an optimum hydrogen flow rate that gives a maximum 
reactor thrust-to-weight ratio. For the conditions of this analysis, an optimum 
will not necessarily occur for all conditions because of the upper limits of 
10 6 pounds for reactor weight and 300 centimeters for cavity radius; in such 
instances the hydrogen flow rate at the cutoff point was used. 
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Turbulent flow . - The effect of hydrogen flow rate on the separation process 
is qualitatively the same for turbulent flow as for laminar flow. Figure 12 
shows the uranium-to-hydrogen atom-density ratio for three hydrogen flow rates; 
separation of the -uranium occurs for the smaller flow rates. 

There is, however, a dissimilarity between laminar and turbulent flow. Al¬ 
though separation does occur for turbulent flow, the actual distribution of ura¬ 
nium in the vortex is such that there is no optimum hydrogen flow rate as shown 
by the radial uranium distributions in figure 13. In the presence of turbulence, 
uranium concentrations are much higher near the outer portion of the vortex than 
for laminar flow. This increase in uranium concentration is a result of the in¬ 
crease of concentration diffusion due to turbulence. Because of this increased 
concentration and the effect of the radial pressure gradients, a decrease of hy¬ 
drogen flow rate gives more separation but actually results in lower average fuel 
density in the reactor. Thus, the "best" performance for turbulent flow is ob¬ 
tained by increasing the propellant flow rate; however, the uranium distribution 
at high flow rates is such that most of the heat generation, and therefore the 
high-temperature region, would occur too near the cavity wall. So, the selection 
of maximum performance for turbulent flow is arbitrary in the sense that it is 
for a hydrogen flow rate that gives a uranium distribution similar to the curve 
in figure 13 for 0.0058 pound per second per foot rather than the one for 
0.01 pound per second per foot. 


Effect of Dissociation 

The double maximums exhibited by the uranium-concentration profiles for tur¬ 
bulent flow shown in figure 13 are an interesting result and deserving of com¬ 
ment. The initial increase of uranium concentration near the outer wall indi¬ 
cates that the separation process has started. The decrease in uranium concen¬ 
tration to a minimum at a radius ratio of approximately 0.65 seems to represent a 
temporary reversal of separation. 

The curve for a hydrogen flow rate of 0.0058 pound per second per foot is 
shown in figure 14 along with hydrogen-species-concentration profiles. The in¬ 
crease in uranium concentration due to the separation process is reversed by the 
dissociation of molecular hydrogen into atomic hydrogen. After the dissociation 
process is essentially complete, the separation process again causes an increase 
in uranium concentration to a second maximum. A s imi lar effect is present for 
laminar flow but is not significant, because the uranium concentration is much 
lower where dissociation occurs than for turbulent flow. 


Illustrative Laminar Example 

The radial variations of the various parameters for a laminar example are 
illustrated in figure 15. Figure 15(a) shows hydrogen molecular weight, tangen¬ 
tial velocity, and static temperature as functions of vortex radius. The static- 
temperature profile is the same as was shown in figure 9(a), and for laminar flow 
the tangential velocity varies inversely with the radius. The hydrogen molecular 
weight is initially at 2 and changes to 1 as dissociation takes place and de¬ 
creases toward l/2 as ionization begins. Other hydrogen properties, such as spe- 
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cific heat, are also taken into account in the program, and these properties are 
discussed m detail in appendix C. 


™ - Plfiure 15 ^i shovs the radial velocity variation of hydrogen and uranium. 

The decrease of the uranium radial velocity at a radius ratio of 0. 6 is a measure 
of the degree of separation; it is essentially the inverse of concentration. 

Hydrogen- and uranium-concentration profiles and the static-pressure varia- 
tion are shown in figure 15(c). The interesting point of this figure is the 
s atic-pressure curve. The decrease of static pressure is moderate until the 
uranium-concentration buildup occurs. The primary decrease in static pressure 
is due to the presence of the uranium cloud. This is typical of laminar flow 

I ne I?^ eparati ° n iS €ncountered - An excessive static-pressure drop is necessary 
to diffuse the hydrogen through the fuel-rich region. 


ROCKET PERFORMANCE CHARACTERISTICS 
Effect of Reactor Length-to-Diameter Ratio 

As was shown in figure 8 , larger reactor cavity length-to-diameter ratios 
give higher thrust-to-yeight ratios. There is some upper limit at which the as¬ 
sumption of axial independence of the flow field would no longer be valid. It 
is not necessary to pursue this point, however, because the gain in performance 
is slight after a length-to-diameter ratio of about 2 is reached as shown in 
figure 16 for the illustrative laminar example. 

The thrust-to-weight ratio goes to 0 at a length-to-diameter ratio equal to 
zero, because there is no hydrogen flow, hut the end moderator-reflector region 
still gives some reactor weight. At length-to-diameter ratios greater than 2 
the increase in critical size causes the reactor weight to increase in proportion 
to the length. 


Effect of Hydrogen-to-Fuel Flow Ratio 

A parameter of economic interest, the ratio of hydrogen-to-fuel flow rates 
affects the performance through the average fuel density (fig. 17 ). The decrease 
in average fuel density with increasing flow-rate ratio is more severe for lam¬ 
inar flow than for turbulent flow. Even for laminar flow, however, this effect 
is less pronounced than the effect that might be anticipated. A flow-rate ratio 
of 100 .to 1 (the ratio used in the performance optimization) gives average fuel 
densities that are only 10 percent higher than those that would be obtained for 
a ratio of 10 , 000 to 1 . 


MAXIMUM PERFORMANCE 
Fuel Distribution 

All the pertinent parameters affecting the separation process were varied 
to determine their effect on reactor thrust-to-weight ratio. The primary effects 
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have been shown and discussed. Tangential velocity was limited to 12, 600 feet 
per second at the outer radius; this velocity corresponds to a Mach number of 
1.0. Vortex-exit radius ratios and hydrogen flow rates were restricted to com¬ 
binations within choking limits. The radius ratio of the maximum mole fraction 
of uranium was varied with the restriction that the uranium-to-hydrogen radial- 
velocity ratio at the outer radius be a reasonable (>10“ 2 ) number. This restric¬ 
tion did not actually limit the optimization, since the best performances were 
for conditions that gave radial-velocity ratios greater than 0.1. 

Figure 18 shows the uranium-concentration profiles for the best laminar and 
turbulent conditions. The performance of turbulent flow could have been improved 
by about a factor of 2 by increasing the propellant flow rate, but this would 
have moved the fuel undesirably close to the reactor wall. The curves of fig¬ 
ure 18 were obtained for a temperature profile with a 15,000° R maximum value. 


Performance Capability 

The parameters for the best performance for both laminar and turbulent flows 
are listed in the following table. For both situations the hydrogen-to-fuel 


Parameter 

Flow 


Laminar 

Turbulent 

Hydrogen-to-fuel flow-rate 
ratio, V]_/wg 

100 

100 

Average uranium density, 
( B 2) av ; particles/cc 

0.46X10 18 

0.55X10 18 

Vortex outer radius, r Q , ft 

5.5 

4.8 

Reactor-cavity length-to- 
diameter ratio, L/D 

2 

2 

Specific impulse, I S p, sec 

3000 

3000 

Thrust, F, lb 

530 

330 

Reactor thrust-to-weight 
ratio, F/V r 

1. 2 3x10“ 3 

<1X10" 3 


flow-rate ratio is 100, the reactor cavity length-to-diameter ratio is 2, and 
the specific impulse is 3000 seconds. The thrust for laminar flow is 530 pounds, 
and that for turbulent flow is 330 pounds. The maximum thrust-to-weight ratio 
for laminar flow is 1.2x10“ 3 and for turbulent flow is less than 1x10" 3 , but just 
how much less is somewhat arbitrary because it depends on how much uranium can be 
maintained near the cavity wall. Additional study of this problem is hardly jus¬ 
tified, however, since the greatest possible improvement for turbulent flow could 
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only yield the performance of laminar flow. 

To reiterate the conclusion drawn from figure 2, a 3000-second-impulse 
rocket must have a powerplant thrust-to-weight ratio of at least 10“ 2 to equal 
the performance of an advanced solid-core system. This ratio is not possible 
with a single-cavity vortex that requires all of the hydrogen to diffuse through 
a fuel-rich region. Since all of the assumptions made in the calculation of the 
performance shown in the table were optimistic, any further effort to determine 
performance capability more realistically must be considered unnecessary. 


ALTERNATE VORTEX SYSTEMS 

Any alternate application of a vortex containment scheme to gas-core reac¬ 
tors must provide an increase in thrust of at least two orders of magnitude by 
either (l) effectively increasing the "length" of the reactor by a multiple tube 
or a. "matrix" arrangement, or (2) flowing only a small (<1 percent) portion of 
the hydrogen through the fuel region and bypassing the remainder through the cav¬ 
ity in such a way that the energy is acquired by thermal radiation. An analysis 
of a multiple tube system is given in reference 1, and a vortex matrix is dis¬ 
cussed in reference 11. 

An illustrative example of a hydrogen, bypass flow configuration is shown in 
figure 19. It should be emphasized that the concept shown is only used to il¬ 
lustrate the principle of bypass flow and does not necessarily represent a feas¬ 
ible arrangement. 

A brief calculation was made to determine the possible gain of such a sys¬ 
tem. The uranium concentration was calculated with the hydrodynamic program for 
the conditions shown. The hydrogen flow rate necessary to achieve this separa¬ 
tion is considerably less than that required to choke the vortex exit. The sup¬ 
plementary hydrogen flow rate was computed from the average vortex-core condi¬ 
tions from the separation program and the choking-limit equation. If the 
tangential-velocity and static-pressure profiles of the supplementary and primary 
vortex flows are matched, the two flows may be superimposed with no effect on the 
vortex strength in the annular region. 

If exit choking is the only limit imposed on the supplementary flow, a reac¬ 
tor thrust-to-weight ratio in the range- of 1 to 10 appears possible. The heat- 
transfer mechanism is now primarily thermal radiation rather than molecular col¬ 
lisions. Because the heat source is between the cavity wall and the bypass hy¬ 
drogen, both the primary hydrogen (to protect the wall) and the supplementary 
flow would have to be seeded in order to absorb radiant energy. The limit on 
performance of such a system will be imposed by the amount of radiation reaching 
the wall and not the separation process. It is possible that such a limit would 
result in only a marginal increase in performance over that of a reactor with 
only primary flow. 

Again, the system shown in figure 19 is only intended to illustrate the 
principle of bypass flow as a possible way to increase the performance capabil¬ 
ities of a single-cavity vortex system. Hydrodynamically, the system appears 
feasible, and the nucleonic characteristics are the same as those investigated 


19 



herein. A precise definition of the performance potential of the general system 
illustrated in figure 19 -would require a detailed radiation heat-transfer anal¬ 
ysis. 


SUMMARY OF RESULTS 


An analysis of a single-cavity vortex gaseous reactor was made for dissoci¬ 
ating ionizing hydrogen gas and uranium 235 nuclear fuel. To obtain the maximum 
performance, as measured by reactor thrust-to-weight ratio, optimistic assump¬ 
tions were made; for example, criticality requirements were used for a reactor 
cavity with an external moderator of Dg 0 100 centimeters thick. An equilibrium 

specific impulse of 3000 seconds was obtained by (l) the assumption that 10 per¬ 
cent of the total heat generation occurs in the solid region and (2) that hydro¬ 
gen enters the reactor cavity at 5000° R. The results of this study are summa¬ 
rized as follows: 

(1) For laminar flow, the maximum attainable thrust-to-re actor-weight ratio 
is 1.2x10-3, hut for turbulent flow a ratio of less than 10" 3 is possible. 

(2) The low thrust results from the condition that all the hydrogen must 
diffuse through the fuel-rich region and that a low flow rate is necessary to 
achieve satisfactory fuel distributions. 

(3) For a typical interplanetary mission, a specific impulse of .3000 seconds 
would require a powerplant thrust-to-weight ratio of 10“ 2 in order to equal the 
performance of a solid-core nuclear rocket operating at an impulse of 1000 sec¬ 
onds. 


(4) For laminar flow, an increase in hydrogen flow rate causes a decrease 
in average fuel concentration in the reactor; consequently, there is an optimum 
hydrogen flow rate that gives a maximum thrust-to-weight ratio. 

(5) For turbulent flow, an increase of hydrogen flow results in an increase 
of average fuel concentration because the fuel is moved nearer to the reactor 
walls; hence, an arbitrary limit on hydrogen flow rate results from wall-heating 
considerations. 

(6) For turbulent flow, the dissociation of hydrogen tends to inhibit the 
increase of fuel concentration due to the diffusion separation process. 

(7) A brief examination of the performance capability of a single-cavity 
vortex system employing supplementary hydrogen flow through the vortex core in¬ 
dicates that a reactor thrust-to-weight ratio in the range of 1 to 10 may be pos¬ 
sible. 


CONCLUSIONS 


(l) The performance potential of a single-cavity vortex gaseous reactor in 
which all the propellant must diffuse through the fuel-rich region is less than 
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that of a solid-core nuclear-rocket engine. 

(2) An alternate scheme that offers sufficiently high thrust-to-weight ra¬ 
tios is one which effectively increases the vortex length and, therefore, the 
propellant flow rate, that is, a vortex "tube bundle" or "matrix" system of many 
large length-to-diameter ratio parallel vortices. Such systems, however, may 
incur severe wall-cooling problems that result from greatly reduced thermal ra¬ 
diation path lengths; if so, this would result in a decrease of attainable spe¬ 
cific impulse. 

(3) A possible application of single-cavity vortex containment to gaseous 
nuclear rockets is to a system wherein a small fraction, for instance 1 percent, 
of the total propellant flow passes through the fuel region. The remaining 99 
percent of the propellant must then bypass the fuel region and acquire energy by 
thermal radiation. Such a system is capable of reactor thrust-to-weight ratios 
of the order of 1 to 10, if it is only limited by the separation process and re¬ 
actor criticality. 


CONCLUDING REMARKS 


The primary conclusion reached herein that a single-cavity vortex gaseous 
reactor where all of the hydrogen must diffuse through the fuel region is inap¬ 
plicable to rocket propulsion, while perhaps neither surprising nor unique, is 
based on a consideration of several effects heretofore treated separately. It is 
expressed or implied, perhaps less specifically, in references 1 and 17, for ex¬ 
ample. Some previous studies of vortex-reactor schemes have included turbulence 
effects; others have considered two-component separation in the presence of heat 
generation, and gas-core nucleonics have also been treated. 

A discussion of the effect of hydrogen dissociation and ionization is pre¬ 
sented in reference 17 but is not included in the separation analysis. To this 
author * s knowledge, the effect of hydrogen dissociation and ionization on the 
separation process has not been quantitatively investigated elsewhere. It Is in¬ 
cluded herein with the optimistic assumption that the uranium does not ionize. 

To this extent, the conclusion presented herein, while unchanged, is on a more 
comprehensive basis. 

The treatment of vortex turbulence in this study follows the general ap¬ 
proach of previous work reported in references 2, 3, 7, 8, 9, 10, and 11. Here 
the deviations of tangential-velocity and static-pressure profiles from laminar 
theory predictions are attributed to the presence of an axially independent tur¬ 
bulent flow field. This procedure appears valid, at least to the extent that It 
has been used to correlate existing vortex data for single-component compressible 
flow. Recent studies (e.g., ref. 20) indicate that at least part of the Indi¬ 
cated deviation from two-dimensional laminar flow results from boundary-layer ef¬ 
fects on the end surfaces of a confined vortex. Such an approach suggests that 
vortex turbulence levels may be less than anticipated. Since a quantitatively 
usable understanding of all of the factors governing single-component vortex flow 
Is not presently available, a more elaborate treatment of a two-component vortex 
than that used In this report is not warranted at this t ime . 
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In order to obtain a maximum, performance capability, it was not necessary to 
consider In detail the effect of uranium-to-hydrogen radial-velocity ratio at the 
point of Injection, (ug/ui) , on ‘tbe separation process. The best performance 

was determined by a free choice of the radius ratio at which the fuel concentra¬ 
tion is a maximum. This choice then determines a value of (ug/ui ) q - a value 

that might, not be physically attainable. So, while the maximum performance is 
readily calculable, a determination of the most probable performance would re¬ 
quire an additional consideration of what boundary values of (ug/u^)^ are pos¬ 
sible in a two-component vortex system. 

A few comments about a bypass flow vortex system are appropriate. Consider 
a single-cavity vortex with a sufficiently low flow rate of primary hydrogen dif¬ 
fusing through the fuel region that a critical concentration exists. If a small 
(i.e., 10 percent of the primary flow) amount of supplementary hydrogen, seeded 
so as to be opaque, Is introduced into the core of the vortex, it will be heated 
by thermal radiation to a temperature that is effectively that of the exhausting 
primary flow. Thus some improvement in performance seems possible. 

In order to obtain thrust-to-weight levels of interest, however, a supple¬ 
mentary flow at least ten times that of the primary flow would be required. The 
fuel-region temperature would have to be increased in order to radiatively heat 
the supplementary flow to a sufficiently high temperature, which, in turn, would 
increase the radiation heat transfer to the reactor wall. It therefore appears 
that an optimum supplementary flow rate would exist In terms of an increased 
thrust-to-weight ratio and a decreased specific impulse. The maximum possible 
thrust-to-weight ratio (l to 10) of a supplementary flow system is high enough to 
warrant a more detailed analysis of the limits imposed by thermal radiation. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 16, 1962 
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APPEJTOIX A 


A 


a l> 

B 

b 

C 

C P 

D 

D 12 

F 


g 

X< 

J 

L 

M 


n 

P 

P 

R 

Re 


r 

T 

t 


u 


SIMBOLS 

area, sq ft 

a 5 numerical constants for temperature profile 
density, particles/cc 
constant in diffusion equation 
parameter in velocity profile 
specific heat, Btu/(lb)(°F) 
reactor-cavity diameter, ft 
binary diffusion coefficient, sq ft/sec 
thrust, lb force 
gravitational constant, ft/sec 2 
specific impulse, sec 

mechanical equivalent of heat, ft-lb/Btu 

reactor-cavity length, ft 

molecular weight, lb/(lb)(mole) 

molar density, (lb)(mole)/cu ft 

total pressure, lb/sq ft 

static pressure, lb/sq ft 

universal gas constant, ft-lb/(lb)(mole)(°R) 

radial Reynolds number 

radius, ft, except as indicated 

total temperature, °R 

static temperature, °R 

radial velocity, ft/sec 

molecular volume of species i, cc/(g)(mole) 
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V 


tangential velocity, ft/sec 
W weight, rb 

w mass flow rate per reactor length, lh/(sec)(ft) 
x vortex radius ratio, r/r Q 

y mole fraction 

z axial coordinate, ft 

a.? parameter in diffusion equation, sq ft 

0 parameter in diffusion equation, sq ft 

T ratio of specific heats 

e eddy diffusivity, sq ft/sec 

0 angular coordinate, radians 

\i viscosity, lb/(ft)(sec) 

p density, lh/cu ft 

cp percent of total power 

^ function in velocity profile (see eq. (ll)) 

Subscripts: 

av average 

c critical 

ch choking 

l laminar 

ma mass 

mo momentum 

n vortex inner radius 

o vortex outer radius 

R reactor 

s sonic 
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t turbulent 

1 hydrogen 

2 uranium 
Superscripts: 

r normalized to value at t 

* indicates radial position of maximum y 2 
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APPENDIX B 


COMPUTER PROGRAM 


By Muriel B. Elan 

At the end of this appendix is a listing of the vortex gaseous-reactor hy¬ 
drodynamic code -written in Fortran for the IBM 7090 digital computer operating 
with the Lewis monitor system. Basically, it is a solution of two first-order 
nonlinear equations (diffusion and momentum for a two-component vortex field) 
using the Runge-Kutta numerical technique to compute uranium mole fraction and 
static pressure as a function of radial position. A block diagram of the pro¬ 
gram is given in figure 20. 

The code can be run for either laminar or turbulent flow ("computed go to" 
controlled). Computation is initiated at an assigned (input) startin g radius 
ratio x*, where dyg/dx is set equal to zero and pressure p is set equal to 

some starting value PRE. The program calculates hydrodynamic parameters from 
x* to x s= 1. The calculation returns to x* and then proceeds to x n . At 
this point, vortex-core equations are substituted for annular equations, and the 
calculation proceeds inward to some arbitrary terminal radius x Q . 

Calculations can be made for either adiabatic or heat-generation (sense- 
switch controlled) conditions. 

When heat is generated, subroutine curve fit (SUB CF) is called. This sub¬ 
routine is based on an arbitrary temperature-distribution curve of temperature t 
against radius x. At each value of x, the static temperature is determined 
from a five-coefficient polynomial in the form of t = + a 2 x + a 3 x 2 + a^x 3 + a 5 x 4 , 

where the coefficients are input to the program in the following form: 


( a l> 

a 2’ 

a 3> 

a 4> 

a 5^]_ 

1 > x X). 725 

( a l> 

a 2> 

a 3* 

a 4, 

a 5) 2 

0.725 > x >0.65 

( a l> 

a 2) 

a 3> 

a 4, 

a s) 3 

0. 65 > x > 0.55 

( a l. 

a 2> 

a 3-> 

a 4> 

a 5> 4 

0.55 > x X) 


Choices of the previous Ax*s are arbitrary and can be varied to fit the curve. 

At each value of x, subroutines SUB INT and SUB YH are called. In SUB INT 
static pressure and static temperature are used to determine hydrogen properties 
(described in appendix C). In SUB YH the hydrogen properties are used to calcu¬ 
late hydrogen molecular weight, specific-heat ratio, mole fraction of H, mole 
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fraction of Hg, and mole fraction of H*. 


When heat is not generated, static temperature, static-temperature ratio, 
otal temperature, and total-temperature ratio are computed from adiabatic-flow 
equations using input values for specific heat, specific-heat ratio, and hydro- 
gen molecular weight. 

An initial increment size of Ax = 0.0005 is used to start the calculation. 
At each subsequent step, the increment size is adjusted according to the test: 


1 'W 4 *! 

Ax 

0-1 

0.0005 

1-10 

0.0001 

10 - 00 

0.00005 


The output is keyed to a counter and increment size so that it occurs at 
Ax intervals of 0.05. 

Accuracy of this problem that uses the Bunge-Kutta numerical technique is 
dependent on the size of the increment. Sizes of increments were decreased and 
problems were run until no significant difference in y 2 was noted. 

For a typical laminar example, doubling Ax resulted in less than 0.1 of 
1 percent change in y 2 near the peak. 


IHPUT 

The program input is required in the following order: 

(1) *Bata 

(2) Binary deck. Table of hydrogen properties (CP, ALPHA, BETA), at 40 

values of temperature for each of 10 pressure levels. (This deck is 
read in once with main deck.) 

(5) Card number 1 Format (14) 

Control words - The "computed go to" words necessary to specify which 
case is to be run are: (l) laminar, (2) turbulent. 

(4) Card number 2 

Sense switch control words are used as follows: Sense switch (l) con¬ 
trols adiabatic (on) or heat generation (off). Sense switch (2) con¬ 
trols computation of eddy diffusivity. Static pressure and temperature 
ratios are set to 1.0 (on). 

(5) Card number 3 Format (7F10.5) 
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W1L 
W2L 
XSTAR 
REMO 
' AM2 
VO 
TO 


hydrogen flow rate, lb/(sec)(ft) 
uranium flow rate, lb/(sec)(ft) 
starting radius ratio (dy 2 /dx = 0) 

Reynolds number for momentum transfer 

uranium molecular weight 

tangential velocity at x = 1.0, ft/sec 

static temperature at x = 1. 0 (for no heat generation) 


(6) Card number 4 Format (7F10.5) 


PRE 

XR 

REMA. 

VIR 

S 

Rim 

CPA 


starting pressure, atm 
vortex-exit radius ratio 
Reynolds number for mass transfer 
molecular volume term in Gilliland equation 
constant for Rewton test 
run number 

specific heat (for no heat generation) 

(7) Card number 5 Format (7F10.5) 

AMI hydrogen molecular weight (for no heat generation) 
GAMMA specific-heat ratio (for no heat generation) 

(8) Cards 6 to 10 Format (4E15.8) 

Constants for curve fit routine 

(9) Card number 11 Format (7F10.5,315) 

^increment step size 

terminal radius 
• constants for output counter 


PELT l] 
PELT 2 M 
PELT 3 J 


XS 


R1 

R2 

R3 



WIL hydrogen flow rate, lb/(sec)(ft) 

W2L uranium flow rate, lh/(sec)(ft) 

XSTAR starting radius ratio (dy 2 /dx = 0) 

REMO Reynolds number for momentum transfer 

AM2 uranium molecular weight 

VO tangential velocity at x = 1.0, ft/sec 

TO static temperature at x = 1. 0 (for no heat generation) 

(6) Card number 4 Format (7F10.5) 

PRE starting pressure, atm 

XN vortex-exit radius ratio 

REMA. Reynolds number for mass transfer 

VXN molecular volume term in.Gilliland equation 

S constant for Newton test 

RUN run number 

CPA specific heat (for no heat generation) 

(7) Card number 5 Format (7F10.5) 

AMI hydrogen molecular weight (for no heat generation) 

GAMMA specific-heat ratio (for no heat generation) 

(8) Cards 6 to 10 Format (4E15.8) 

Constants for curve fit routine 

(9) Card number 11 Format (7F10.5, 315) 

PELT 1*1 

PELT 2 > Increment step size 
PELT 3 J 

XS terminal radius 

HI *1 

N2 > constants for output counter 
N3 J 
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OUTPUT 


Input parameters, type of calculation, and starting values at x* are the 
initial output. Then pertinent parameters are read out at intervals of approxi¬ 
mately 0.05 in x. 


RUNNING TIME 

Program running time is somewhat dependent on the input, since Ax is ad¬ 
justed as a function of dyg/dx; for heat-generation runs the average running 

time on the IBM 7090 with the Lewis monitor system is approximately 3 minutes. 
The no-heat-generation runs, because hydrogen-property interpretations are not 
required, have a running time of about 2 minutes. 
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VORTEX GASEOUS REACTOR HYDRODYNAMIC PROGRAM 


DIMENSION T(10*40)*CP(10*40)*ALPHA(10*40)*BETA!10*40)*P(10) 
DIMENSION AT(40)*ACP(40)*ALP{40)*BET < 40)*PPP(4)*TEMPI(4) 

DIMENSION TEMP2(4)*TEMP3(4) 

COMMON T*CP*ALPHA * BETA*P,AT*ACP*ALP*BET»PPP*TEMP1*TEMP2*TEMP3*PP*T 
IT *ANS1*ANS2*ANS3*YH2 *YH*YHPLUS*GAMMA*AMI*X* ZA.ZB*ZC*ZD»ZE*YA,YB* 
2YC*YD*YE*XA*XB*XC*XD*XE*WA*WB*WC*WD*WE*XSQ*X3 *X4 *TLIL 
CALL BC READ(T(10*40),T( 1 * 1 )) 

CALL BC READ(CP(10*40)* CP(1 * 1) ) 

CALL BC READ(ALPHA(10*40)*ALPHA{1*1)) 

CALL BC READ(BETA(10*40)*BETA(1*1)> 

CALL BC READ{P(10)*P{1)) 

3 FORMAT(4F10.5*315) 

4 FORMAT{1H0E15•8 *6E17* 8) 

6 FORMAT(1H0*9X * 2HCP*13X * 5HALPHA»11X.4HBETA *12X*4HTEMP »13X*4HPR£S) 

7 FORMAT(E15*8 *6E17*8 ) 

8 FORMAT(E15*8*E17«8*E16#8*E17*8*E16*8*E17«8#E16*8*F6*0) 

9 FORMAT(7F10#5) 

10 FORMAT(4E15*8) 

33 FORMAT(21H0 PRESSURE OVER 1000*) 

39 FORMAT(1H0*7X*2HM1*13X*5HGAMMA*12X*2HYH*14X*3HYH2*13X•6HYHPLUS*13X 
1 * 2HCP) 

40 FORMAT(1H0E25.8*E29.8 *E28.8 *E16.8 *E16.8) 

41 FORMAT(106H0 X Y2 

1 U2/U1 B2 C2) 

42 FORMAT(1H0*8X *1HV*12X *6HTS/TS**12X*2HB1*14X*2HC1*12X*6HTT/TT**10X» 
?6HPS/PS**10X*6HPT/PT**11X) 

43 FORMAT(20H0 NO HEAT GENERATION) 

44 FORMAT(17H0 HEAT GENERATION) 

45 FORMAT(16H0TLIL MINUS STOP) 

46 FORMAT{16H0U2U1 MINUS STOP) 

47 FORMAT(16H0 Y2 MINUS STOP) 

48 FORMAT(18H0 Y2 OVER ONE STOP) 

49 FORMAT(E15*8*7E16*8) 

71 FORMAT(1H0*7X*1HV*13X*6HTS/TS**12X*2HB1 *14X *2HC1*10X*6H TEMP*12X* 
16H PRES) 

95 FORMAT(14) 

97 FORMAT f 62H1 
1R0GRAM) 

98 FORMAT{62H1 
1R0GRAM) 

99 FORMAT{102H0 PEAK—TO—WALL CALCULATIONS FOR 

1VORTEX CAVITY REACTOR RAGSDALE-EI AN PAX 2115) 


TURBULENT P 


LAMINAR P 


106 FORMAT(118HO 
1 REMO 

107 FORMAT(120H0 
1 TLILS 

108 FORMAT(120H0 
1 PS* 

109 FORMAT(62H0 
1 ) 

PI=3.14159 
G*32.2 
AJ=778•0 


W1L 

AMI 

PSI1 

TBIGS 

ALPH2 

REMA 


W2L 
AM2 
PS 12 
BP 
Y2STA 

XN 


XSTAR 

VO 

PSI20 

D 

U2U1S 


RUN) 


RG=1544* 0 

22 READ INPUT TAPE 7*95.IG 
GO TO (208*207,207),IG 
208 WRITE OUTPUT TAPE 6*98 
WRITE OUTPUT TAPE 6*99 
GO TO 2 

207 WRITE OUTPUT TAPE 6*97 
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WRITE OUTPUT TAPE 6*99 
2 CALL READ SS 

READ INPUT TAPE 7 * 9*W1L*W2L>XSTAR#REMO *AM2*VO *TO*PRE*XN*REMA*V IN►$ 
1 *RUN * CPA *AMI *GAMMA 

READ INPUT TAPE 7*lO*ZA*ZB*ZC*ZD*ZE*YA*YB*YC*YD*YE*XA*XB*XC*XD*XE* 
1WA*WB*WC*WD*WE 

READ INPUT TAPE 7 *3*DELT1*DELT2*DELT3*XS*N1*N2*N3 
111 KTR=0 
RE=REMO 

C=(l./(2*~RE))*(RE/(EXPF(RE/2.)-l.)) 

VOROP=l•0+C*(XN**(RE-2*)-*l.) 

TEM=(1•O/VOROP)**2 
TEM1=(1.0-C)**2 
TEM2=C—C*C 
TEM3=(XSTAR/XN) 

T£M4=TEM3**(RE-2#) 

TEM5=TEM3**(2.*RE-4*) 

TEM6=TEM2/TEM4 

TEM7=C*C/TEM5 

TEM88=TEM1+2.*TEM6+TEM7 

PSI1=TEM88*TEM 

V=(VO*SQRTF(PSIl))/XSTAR 

TEM89=TEM1*( 1 *-4* / RE) “*( 4* /RE*TEM6)-TEM7 

PSI2=TEM89*TEM 

XSQST=XST AR*XST AR 

X3ST=XSTAR*XSQST 

X = XST AR 

XSQ-XSQST 

X3=X*XSQ 

X4=XSQ*XS0 

TEM 90=2•*CPA*G*AJ 

TEM91=VO*VO 

IFCSENSW(l)*503*504*504 

503 PSI20=TEM*(TEM1*(l.-4./RE)-*(<4*/RE)*<TEM2/<1*/XN)**(RE-2•Hl-(C*C/ 
1(1./XN)**(2.*RE-4.H ) 

WRITE OUTPUT TAPE 6*43 
TINFI=T0+PSI20*TEM91/TEM90 
TLILS=TINFI-(PSI2*TEM91/(TEM90*XSQST)) 

IF(TLILS)50.51*51 

50 WRITE OUTPUT TAPE 6*45 
GO TO 22 

51 TBIGS=TLILS+V*V/TEM90 
GO TO 513 

504 CALL CF 
TT=TLIL 
PP=PRE 
TLILS=TLIL 
CALL INT 

WRITE OUTPUT TAPE 6*6 

WRITE OUTPUT TAPE 6*7 *ANSI♦ANS2*ANS3*TT*PP 
TEM 9G=2#*ANS1*G*AJ 
TBIGS-TLILS+V*V/TEM90 
WRITE OUTPUT TAPE 6*44 
CALL SUB YH 
513 TEM92=G*RG*TLILS 

BP=fAM2-AM1)*TEM91*PSI1/TEM92 

DN=(2•62E—6*SQRTF{(l./AMl+l./AM2)*TLILS))/VIN 
C0N8=W1L/W2L*(AM2/AM1) 

B=1.0+CON8 
CON=W1L+W2L 
C0N1=2**PI*REMA 



TL2=DN 

ALPH2=W2L/<2#*PI*AM2*TL2) 

CON5=ALPH2*XSQST/BP 
CON6“CON5*B—l•0 

C0N7-SQRTF < CON6*CON6+{4«*CON5 )) 

Y2STA=((-C0N6 )+C0N7)/2*0 

407 FY=(-ALPH2*(B*Y2STA-1.0)/XSTAR > + (BP*IY2STA-Y2STA*Y2STA)/X3ST) 
IF (ABSF (FY ) “*( S#Y2STA1) 405 *405 *406 
406 FPY-(C-ALPH2*B)/XSTAR)+(BP/X3ST)-(t2*0*BP/X3ST)*Y2STA) 
DY2=FY/FPY 
Y2STA=Y2STA-DY2 
GOTO 407 

405 U2U1S=(a*/Y2STA)~l.)/CON8 
POPST = 1•0 
P0P0=P0PST 
X = XST AR 
Y 2=Y2ST A 
Yl=1•0-Y2STA 
TEM15-GAMMA/{GAMMA—l*0 ) 

TEM16=XSQST*GAMMA*TEM92 
TEM17-(GAMMA-1•0)/2•0 

TEM18=(TEM91*PSI1*(Y1*AM1+Y2*AM2>)/TEM16 
P4S=P0P0*(1•0+TEM17*TEM18)**TEM15 
PP=PRE 

WRITE OUTPUT TAPE 6*106 

WRITE OUTPUT TAPE 6 *8 * W1L *W2L * XSTAR *REMO* AMI * AM2 . VO ■♦RUN 
WRITE OUTPUT TAPE 6*109 
WRITE OUTPUT TAPE 6*7*REMA*XN 
WRITE OUTPUT TAPE 6*107 

WRITE OUTPUT TAPE 6*7 * PS 11»PS12 * PS I 20*TLILS*TBIGS*BP*DN 
WRITE OUTPUT TAPE 6*108 

WRITE OUTPUT TAPE 6*7*ALPH2*Y2STA»U2U1S*PRE 
WRITE OUTPUT TAPE 6.41 
IFCSENSW(l))64,63*63 
64 WRITE OUTPUT TAPE 6*42 
GO TO 12 

63 WRITE OUTPUT TAPE 6*71 
WRITE OUTPUT TAPE 6*39 
12 TEM3=X/XN 

TEM4=TEM3**(RE-2#) 

TEM5=TEM3**(2.*RE-4«) 

TEM6=TEM2/TEM4 

TEM7-C*C/TEM5 

TEM88=TEM1+2**TEM6+TEM7 

PSI1=TEM88*TEM 

V=(VO*SQRTF(PSI1)>/X 

TEM89-TEM1* 11•—4*/RE)-(4*/RE*TEM6)-TEM7 

PSI2=TEM89*TEM 

XSQ-X*X 

IF(X-l*0000001)525*525*22 
525 IF(SENSW(l))515*516*516 

515 TLIL=TINFI-(PSI2*TEM91/tTEM90*XSQ)) 

TTLIL=TLIL/TLILS 

TBIG=TLIL+V*V/TEM90 

TTBIG=TBIG/TBIGS 

GO TO 517 

516 X3=XSG*X 
X4=XSQ*XSQ 
CALL CF 
TT^TLIL 
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TEM 90=2**ANS1*G*AJ 
T8IG=TLIL+V*V/TEM90 
TTLIL=TLIL/TL1LS 
TTBIG-TBIG/TBIGS 
IF(TLIL)50*517*517 
517 TEM92=G*RG*TLIL 
PlOPO=Yl*POPO 
P20P0=Y2*P0P0 
IF{SENSW(1H66*67*67 
67 CALL INT 

IF(PP—1000#) 83*83*84 
84 WRITE OUTPUT TAPE 6*33 
GO TO 22 
83 CALL SUB YH 

66 BP=(AM2-AM1)*TEM91*PSI1/TEM92 

DN = ( 2•62E-6*SQRTF((1./AM1+1•/AM2)*TLIL)J/VIN 
C1=2*22E-3*RG*TLIL*W1L /(2.*X*AM1*P10P0>*PRE*14.7 
C2 = 2 ♦22E“3*RG*TLIL*W2L /(2.*X*AM2*P20P0)*PRE*14.7 
Bl=(1.385E+24*P10P0)/(RG*TLIL)*PRE*14.7 
B2=(1.385E+24*P20P0)/(RG*TLIL)*PRE*14#7 
TEM 15=GAMMA/(GAMMA-1*0) 

TEM16=X*X*GAMMA*G*RG*TLIL 
TEM17=(GAMMA-1•0)/2 *0 

TEM18=(TEM91*PSI1*(Y1*AM1+Y2*AM2))/TEM16 
P4=P0P0*(1*0+TEM17*TEM18)**TEM15 
PPBIG=P4/P4S 
C0N8=W1L/W2L*(AM2/AM1) 

U2U1=((1*/Y2)~1*0)/C0N8 
54 GO TO(526 * 527 * 527)* IG 
527 TL2=DN 

IF (SENSW(2))634>635*635 

634 POPST-1*0 
TSTOT=l•0 
CON6=2**PI*REMA 

EN=(CON*POPST*TSTOT)/(CON6*tY2*AM2+Y1*AM1)) 

GO TO 530 

635 CON6=2«*PI*REMA 
C0N9=TLILS/TLIL 

EN=(CON*POPO*CON9)/(CON6*(Y2*AM2+Y1*AM1 )) 

GO TO 530 
526 TL2=DN 

530 ALPH2=W2L/(2.*PI*AM2*TL2> 

GO TO (532>533*533)*IG 

532 TEM8=1•0 
GO TO 52 

533 TEM8=,DN/ (DN+EN ) 

52 DYDX =(((Y2-Y2*Y2)*BP/X**3)-((B*Y2~1*0>*ALPH2/X)>*TEM8 
Y2P=SQRTF(DYDX*DYDX) 

IF(Y2P-1.0)13*13*14 

13 D£LTX=DELT1 
N=N1 

GO TO 17 

14 IF(Y2P-10.)15*15*16 

15 DELTX S DELT2 
N=N2 

GO TO 17 

16 DELTX=DELT3 
N=N3 

17 IF (XMODF(KTR*N))18*19*18 
19 KTR=KTR+1 

WRITE OUTPUT TAPE 6*40*X♦Y2*U2U1*B2*C2 



IF fSENSW(l))80#81#81 

81 WRITE OUTPUT TAPE 6 *7 *V*TTLIL•B1*C1* TT»PP 

WRITE OUTPUT TAPE 6 t4 *AMItGAMMA tYH *YH2 >YHPLUS #ANSI 
GO TO 20 

80 WRITE OUTPUT TAPE 6f49,V*TTLlL*Bl»Cl»TTBIG»PP»PPBIG 
GO TO 20 
18 KTR=KTR+1 
20 V02=TEM91 

TEM13=G*RG*TLIL 

AAK1=TEM8*(—ALPH2*(B*Y2-1*0)/X+BP*(Y2—(Y2*Y2) ) /X**3)*DELTX 
AAL1=(((1*-Y2)*AM1+Y2*AM2)*V02*PSI1*P0P0*DELTX)/(TEM13*X**3) 
TEM9=X+(DELTX/2*0) 

TEM10=TEM9**3 
TEMll=Y2+(AAKl/2*0) 

AAK2=TEM8*((-ALPH2*(B*TEM11—1*0)/TEM91+(BP*t TEM11—TEMI1*TEM11) /TEM 
110 ) )*DELTX 
TEM12=Y2+(AAK2V2»0) 

AAL2=(((1•—TEM11)*AM1+TEM11*AM2)*VQ2*PS11*{POPO+AAL1/2•0)*D£LTX)/( 
1TEM13*TEM10) 

TEM15=X+DELTX 

AAK3=TEM8*<-ALPH2*{B*TEM12“1.0)/TEM9+BP*(TEM12--TEM12*TEM12>/TEM10> 

1*DELTX 

TEM14=Y2 +AAK3 

AAL3= f(l.-TEM12)*AMl+TEM12*AM2)*V02*PSU*(POPO+AAL2/2*0)*DELTX/(TE 
1M13*TEM10) 

AAL4=((l.O-Y2+AAK3)*AMl+AM2*(Y2+AAK3n*V02*PSIl*{POPO+AAL3)*DEL 
1T X/{T E M13 *(X +DELTX)**3) 

AAK4=TEM8*( -ALPH2*(B*TEM14 ~1*0>/TEM15+BP*{TEM14~TEM14*TEM14J/TE 
1M15**3)*DELTX 

DELTY2=(AAKl+2.0*AAK2+2*0*AAK3+AAK4)/6.0 
DP0P0=(AAL1+2#0*AAL2+2.0*AAL3+AAL4)/6*0 
Y2=DELTY2+Y2 
IF(Y2 )55 * 56* 56 

55 Y2-0 

56 IF ( Y2—1*0000001 ) 57* 57 »58 

58 WRITE OUTPUT TAPE 6*48 
GO TO 22 

57 Yl = l# 0-Y2 
P0P0=DP0P0+P0P0 

IF (SENSW 11H 73 *74* 74 

73 PP=P0P0 
GO TO 76 

74 PP=POPO*PRE 
76 X-X+DELTX 

IF(X-l•0000001)59♦5 9♦60 

59 GO TO 12 

60 GO TO 22 
END 
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APPENDIX C 


HYDROGEN PROPERTIES 


Hydrogen.properties used in this study were computed from results presented 
in reference 19. The properties of interest herein are (l) equilibrium species 
mole fraction of H 2 , H, and H 1 "; (2) specific heat at constant pressure; (3) ratio 

of specific heats; and (4) average hydrogen molecular weight. The symbols used 
in this appendix are defined in the following paragraph and, where applicable, 
are the same as in reference 19. 

For a range of temperatures (°K) and pressures (atm), reference 19 gives a 
tabular listing of a, J3, and H/M, where H/M is enthalpy in calories per gram. 
For N initial atoms of H, the equilibrium species concentrations are, by defin¬ 
ition: 

aN/2 Hg molecules 

pN E* protons 

(3N electrons 

(l - a - 0) H atoms 

From these parameters of reference 19, the values used herein were computed 
from the following equations: 

_ _a_ 

yH 2 2 + 20 - a 


y H * 


1 - a - J3 
1 + -P - (a/2) 


P 

1 + 0 - .(a/2) 


M&v = 2 * 016 XH 2 + l-008(y H + y^) 


^(h/m) 

at 






jp *1 - t 2 

where c is at a temperature equal to (t-j_ + tg)/? 


T 


c p ^ JCpMa-y/B 


( Jc p M L/ R ) ' 1 
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The input to the hydrodynamic program properties subroutine was in the form 
of tables There were ten tables, one for each of the following pressures: 1CT 5 
10 , 10 , 10 , 10 , 1, 10, 100, 300, and 500 atmospheres. For each pressure, 

values of a, p, and Cp were listed for 40 temperatures ranging from 900° to 

39, 600° R. A five-point two-variable Lagrange interpolation procedure was used 
to obtain values of a, p, and Cp for any desired pressure and temperature con¬ 
ditions in the vortex. 

The output of this subroutine is shown in figure 21 for pressures of 1, 100 
and 500 atmospheres and for temperatures from 5000° to 22,000° R. Figure 21(b) ’ 
the average hydrogen molecular weight, represents the most significant effect of 
dissociation and ionization of the propellant on the separation. The computed 
average molecular weight affects both the binary diffusion coefficient and the 
static-pressure gradient of the hydrogen in the vortex. Equilibrium specific 
heat is of considerable importance, of course, but for the arbitrary temperature 
profile employed in the hydrodynamic program this effect is hidden. 
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Figure 3. - Effect of specific impulse on mission 
time for seven-man Venus mission. Initial 

•weight in orbit, 1.7X10 6 pounds; thrust per 

reactor weight, >10 
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Figure 4. - Gaseous nuclear rocket specific-impulse - temperature relation 
as a function of reactor pressure and hydrogen enthalpy fraction gained 
in solid region at 5000° R inlet temperature. 
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Lgure 6. - Comparison of critical densities for three types of externally 
moderated cavity reactors. 
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Figure 7. - Reactor (moderator) weight as a function of 
reactor-cavity radius and cavity■length-to-diameter ratio 
Moderator thickness. 100 centimeters. 
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(a) Temperature profiles. 

Figure 9. - Effect of temperature profile for laminar flow 
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Figure 10. - Effect of hydrogen flow rate on separation for laminar flow. Hydrogen- 
to-uranium flow rate, 100; inlet static pressure, 500 atmospheres; inlet total 
temperature, 5000° R; peripheral Mach number, 0.9. 
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Figure 11. - Effect of hydrogen flow rate on re¬ 
actor thrust-to-weight ratio for laminar flow. 
Hydrogen-to-uranium flow rate, 100; peripheral 
Mach number, 0.9; inlet static pressure, 

500 atmospheres; inlet total temperature, 

5000° R; reactor-cavity length-to-diameter 
ratio, 2. 
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Figure 12. - Effect of hydrogen flow rate on separation for 
turbulent flow. Hydrogen-to-uranium flow rate, 100; inlet 
static pressure, 500 atmospheres; inlet total temperature, 
5000° R; pheripheral Mach number, 0.9. 
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Figure 13. - Effect of hydrogen flow rate on uranium distribu¬ 
tion for turbulent flow. 
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Vortex-radius ratio, x 

Figure 14. - Effect of dissociation on turbulent 
separation process. 
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(a) Variation of tangential velocity, molecular weight, and 
static temperature. 


Figure 15. - Illustrative laminar example. 
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Figure 16. - Effect of reactor-cavity length-to-diameter ratio 
on reactor thrust-to--weight ratio for illustrative laminar 
example. 



Figure 17. - Effect of hydrogen-to-uranium flow ratio on average fuel density. 
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Figure 21. - Concluded. Equilibrium properties of hydrogen gas. 
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